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ABSTRACT 


Visibilivy is sone of the most important weather elenencva 
jm marine navigation. It is commonly believed that the pri-= 
mary factors determining visibility are the wind speed and 
the relative moisture content of the air. Using regression 
analysis, an attempt is made to establish a linear functional 
relationship between observed visibility and objectively- 
@ectermined parameters from.conventional marine surface obser= 
menos. forecast visibilities are obtained from derived 
regression equations using elements available from a numeri- 
fe primitive equation forecast model. 

Twenty-eight hundred observations in the North Atlantic 
memeenanalyzed,. It was determined that reliable regression 
eauavions were derived only from those observations obtained 
from ships with trained meteorological observers on board. 
mime results indicate that numerical forecasts of surface 
memperature and humidity have not yet reached the accuracy 
Meeesoary for use in this type of statistical prediction 


scheme. 
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i eR ODUC men 


Visibility is one of the most important Weaver = lena. 
i marine navigation but it is one of the most compi@icated 
ofall meteorological observations. The visual range is 
defined as the distance that an object can be seen fi] ‘ 

The horizontal visual range near the surface is called the 
Power lives Semevailingesvastoility is defined @s the hipit@sr 
visibility that is equalled or exceeded over sectors of the 
weed zZon Which; when combined, total one-half or more of the 
Merizon circle LE Preevalling viszbility is Uhe™parancver 
Mmeported by ships at sea. 

Huesca, the absence of objects at known distances makes 
mee ODServation of visibilitvyehnrehly subjective. There, 
M=maopeamancenOf theshorizon, aseobserved from different 
Merels, yieclas the west estimate of the prevailing viswoidic.,. 
The experienced seaman knows how difficult it is to estimate 
visibility at night. The presence of a "loom" around the 
Ship's navigation lights is frequently an indication of 
meeeriorating visibility. It has been estimated that the 
vVasibility in full moonlight is about 20% of the visibility 
in daylight in the same atmosphere Ele 

There are several possible approaches to the problem of 
analyzing and forecasting visibility at sea on a large scale. 
Presently, visibility forecasts are prepared subjectively 


oy 5 o ‘ “ ; ~ s+ = 4 : 
a=C are basec on ship reports, air mass characves, cer 





present and past weather, and climatology, —Unicriy ee 
results are highly dependent upon the experience Teveluon 

the forecaster. Some progress has been made in the area of 
modeling. However, the physical-dynamical models are based 

On emaicrophysmes and requimee input parameters on a Seale thar 
feemot avallable from conventional marine surface observa— 
peems. in this case, numerical generation of data from exist— 
me Observations could be utilized. 

Another approach seems worth trying for immediate opera- 
peeonal use. A simple method is the application of statista— 
cal regression analysis based on objectively-determined ob- 
served and forecast parameters. This is the approach used 
in this paper. An attempt will be made to establish a linear 
functional relationship between the prevailing visibility 
and objectively-determined parameters available from conven- 
Mieaal marine surface observations. Further, an extensi?on 
Will be made to obtain forecast visibility from elements 
Zomelapile {from numerical or other objective prognostc 


techniques. 





II. BACKGROUND 


A review of the available literature reveals much Worker 
peocememal nature about fog prediction. ~However eta cre 
iiGele research on the subject of sea fog was found and few 
attempts had been made to use a statistical approach. 

Schramm attempted to relate Fleet Numerical Weather Central's 
vapor pressure difference fields (saturated vapor pressure- 
ambient vapor pressure) and other available parameters to 
observed visibilities and it is his report that forms the 
basis for this paper [ 4). 

The North Atlantic area was chosen as the locale to per- 
form the analysis. This choice was based on climatology, 
Memoaty of observations, and the availability of previous 
meoearch for comparison wale June and December 1971 were 
the selected periods of interest as climatology indicated 
eemieh percentage of low visibilities during these months. 
Input data for this research were obtained from the Fleet 
memerical Weather Central climatology files. 

The statistical resources of the W.R. Church Computer 
Center at the Naval Postgraduate School are many and varied. 
A stepwise linear regression program called BMDO2R was selec- 
ted for the analysis [6]. Stepwise regression analysis 
enables the user to "build up" his regression equation from 
a Simple linear to a multiple linear form by introducing one 
independent variable at a time. 


ee 
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iis» PROCEDURE 


A. DATA PROCESSING 

ine area of interest was furtmer subdividedeepac amen 
Mmmatology, for each of the two time periods. Figure I 
mmistraves these areas. The three general areas e€mosen 
Mime tetecrred to nerein as the Grand Banks area yw tne dere 
Beetan sea area, and the Station Ship area. 

As the observations from the Ocean Station Vessels alone 
mel trol COMprise a lame enone data sample for each loca- 
men, all reports within two degrees of the assigned statwenm 
mre DOsIulon were utilized. This was regrettable injvuhay 
the Station Ship observations were desirable as a "control" 
iia, LO be compared with other observations in the same 
Beneral area. 

Hearimessuriace observations for Juneyand December 7 (@ 
Bemenwcxtracted from the Fleet Numerical Weather Central 
memauelopy files. Whese raw data are stored as received 
ena are not processed prior to being stored on magnetic tape. 
femeess error check and rejection of duplicate observavioens 
Mere first performed on the data. Next, the following para-— 
MeLrers were extracted: 

ie) eee diodl hy eves?) 
2. Wind direction (Dir) 
3. Wind speed (V) 


4 A - ic Sree) Saimna ry 
ned e x1Yr Re — apg Sal or gl L we oS 
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Sea surface temperature (T_) 
Dew point temperature (T) 


Present weather code (WW) 


CoO NN OF UW 


. Total cloud covers) 


\O 


Latitude (LAT) 
10. Longitude (LONG) 


11. Date (DAY) 


The following parameters were generated from the basic data 


above: 


1. eastward component of the wind velocity (u) 
2. northward component of the wind velocity (v) 
. dew-point spread (T-T,) 
-. air-sea temperature difference CEE) 


vapor pressure of the air Ce.) 


3 
4 
5 
6. Saturation vapor pressure of the air (omy 
T saturation vapor pressure of the sea surface Ce) 
8 saturation vapor pressure difference (e em 
9 vapor pressure deficit (eRe) 

10. relative humidity (RH) 

lie WevT_ fOrwetatl Cnesitps OnLy (V-VE_) 


12. evaporation (E) 


13. total water vapor (WVT) 
The vapor pressures were generated from Teten's expression 
[7]: 


AT 


e(mb) = 6.11 eto © ky 


aS 





where T is temperature and A and B are constants. Dewpoint 
(T.) and sea surface temperature re) were entered for vapor 
pressure of the air and saturation vapor pressure at the 

sea surface, respectively. Laevastu [8] developed the fol- 


lowing Dalton-type expression for evaporation: 
ie (e,, - e,? is positive: 


Bin) = (0.20 +ecron;, Vv) Ce, ~ e,) (2) 


Ibe oe - e,) is negative: 


E(m) = (0.077V) (é, - e,) (3) 


a e 


Total water vapor [9] was computed from a regression equa- 
tion relating surface dew-point to total water vapor, shown 


below: 


In WT = -1.288 + 0.03847, | (4) 


The Fleet Numerical Weather Central sea surface temperature 
emalyses* were wsedmas imput to obtain*walwes of the gradient 
©m@ the sea surface temperatures. Values of this gradient 
mae then interpolated to the station ship positions and 
Wovt. Computce. 

Mil parameberseexcept pressure terms were convervedure 
the meter-kilogram-second (MKS) system of units. Over 2800 


observations were utilized. 


B. ANALYSIS OF THE DATA 


Regressivon analysis is one of the most pepularstacis— 


4 


moaleeeols in use. 


U/ 


tepvise resression eCnae ten] 
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mesearcher to vexpleaperie power and speed of the ceG@impmeer 
to examine his data in much greater detail than senoaea. 
multiple linear regression. 

BMDO2R computes a sequence of multiple linear regression 
equations ina stepwise manner. At each step one variable 
jmomagded £O the regression equation. The variable added 
men Lhe sone whach makes the greatest reduction any neem 
ei vie errors squared. At each step of the preeess cempteve 
aaLormation is presented about the current regression equa-— 
tion. Initially, the program accepts up to 99 independent 
variables (predictors) and up to 9,999 cases. 

Bor this study, the dependent variable, visibility, 


was entered in three forms: 


1. Veisibility in meters, 
Co Sibel ty coded stiremeship weportsmeand 


3. £Neturaik log GF unewvisibi lity code. 


The independent variables were separated into two 


groups: 


1. <All available observed parameters (except pre- 
sent weather and cloud cover), and 
2. parameters forecast by the Fleet Numerical 


Weather Central primitive equation model. 


~ 


However, upon examination of the data from the primivave 
equation model, it was found that the forecasts of air tem-— 
perature and vapor pressure were not useable, especially 


Cine  Cerlods Of thern 


mM 


Sane te he a , se 4 
1 instability. ~Censequen as 


IS 





analyzed fields of surface temperature and vapor pressure 
were obtained. 

No attempt was made to separate the observations with 
precipitation causing reduced visibility. Had this been 
done, the reduced sample size would have precluded any sig- 
“nificant results. 

A regression analysis was run for each area and for 
all combinations of observed see ecenaeee variables and the 
three forms of the visibility. Forecast visibilities were 
mmern OOLained by uSing thé regression equations derived from 
the observed parameters. Two "forecast" visibilities were 
generated, one without the unuseable forecasts of air temper- 
ature and vapor pressure, but including all other parameters. 
The other "forecast" visibility entailed the use of analyzed 
fields of air temperature and vapor pressure with the re- 


maining forecast parameters. 
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IV. “RESULTS 


Tables |] through 11 iliustrate the Best  ,eSmiirs serene 
meune observed parameters for each area during the two 
fme periods. The total number of observations used dn 
the analysis for each particular area is shown. The coef- 
ficient of determinancy, R*, can also be considered as the 
mreactCion of variance explained. In all cases an b=ratio 
Paeaver thanwc.25 indicates a level of Sipnaficance preaver 
than 0.95. Note that the standard error of estimate is 
mupmessco wcliner in meters or units of the visdbi lary ween 
time most significant variables are those that caused the 
greatest reduction in the sum of squares of the errors 
moumare liseed in decreasing order of importance. The 
@erived r€gression equation with coefficients is also 
listed. The parameters are expressed in the MKS (meter- 
kilogram-seconds) system of units except for temperature 
(°C) and pressure (millibars). 

Dmecencmad, less thammone-half of the variance between 
mae observed visibility and the visibility computed by 
the regression equation is explained. However, in nearly 
every case the F-ratio is large enough to indicate that 
a significant relationship exists. As an average of all 
eases, tUne=noltowing input variables are 1listeqmin oman 
Of importance as the most significant in reducing enews 


re Fa ~ 
PeooWeteo Of ule "CrrcCran 
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TABLE I 
Area of Station Ship A 


62N 33W 
a ET 
=a ae 
a O S76 0.5174 


Most Significant 


él Tl Ria—T.,e RH, U—T. , Vice ge 
Variables : q's d 





Regression Equations 
jume: Vis = 130m068 - Of leten — C20292c7 
~41320T - 0.049457. + 2.16727e, = SOB Se 
-07973e,, = 0. G412ISE = 0 737e52nes 1551393 VevT. 
.0O1211Dir - 0.08479V - 0.48006(T-T,). 


December: Vis = 98.864 - 0.07318u + 0.O0H4174v 


uso Rou — 0.62250T. - 2.5038e, + Blea S 


sie 262he— >. 313N- V0 = OC l222ber — Oe e 
,66596(1-T,)- 
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TABLE II 
Area of Station Sshaipee 
56.5N 51W 


Number of 70 35 
. Observations ; 


— 


Most Significant v,e,-e.,T-T E,V,T-T, 
Variables 


Regression Equations 
June: Vis = —- 624933 - 74.8296u - 375.770v 
+6687 .87T. i Bie Sf See) - 10482 .2(e,-e_) - 485.580E 
+ 7003.66RH - 171.426Dir - 203.656V + 38632.8(T-T.) 


December: Vis = 95.8679 + 0.02769u + 0.01104v 
- 0.499357 + TEESE. ~ O28 1 Gewae ~) + 0.60993E 
- 0.14112RH + 45.8342wWVT - 0.01723Dir - 0.21261V 


- 0.39344(T-T.). 





Ie 





WABIGE, Jie 
Area of Station Ship C 


Dea (Ne so om 


RES = December 


Number ene 76 63 
Observations 


Regression Equations 


June: Vis = 143.975 - 0.62642u - 0.19038v 


1.04475T + see ae Te 2 eB es - 15.3000e,, 


Oo Oo Cenae. ) + 0.81175E - 0.07947RH + 148.028V-VT. ; 


0.19316Dir - 1.28930T, 


December: Vis = - 54141.4 + 27.7416u + 62.3268v 
2251.00T. - 2433.96(7_-T) ~- 14665.7e, + 7 1GGR23 (e Saem) 
1784.28E - 516.637RH + 769055.WVT + 82.1492Dir 
§9.6531V. 
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TABLE IV 
Area of Station Ship D 
HUN 41W 


° 


Number of 81 101 
Observations 
O< 7000 0.6747 
0.4901 0.4552 


F-Ratio 5.446 
Standard Error 4777 meters 5057 meters 


Mest Significant 
Variables 












RH, Tae 


Regression Equation 


June: Vis = 159594. - 60.0878u - 156.446v 
+ 942.632T_ - 6494.18(T_-T) + 2891.01e, - 2848.22e,, 
+ 7416 .06(e,-e.) + 306.969EF - 1672.83RH + 27.3565Dir 
= 245.243V - 4486.08(T-T,) 


December: Vis = 45958.3 - 50.2461u + 114.716v 
- 6863.42T. - SDE ALOE Se) - 3071.13e, + 10147 .8e_, 
- SCC eee cs) —~ 228.604EF - 541.377RH + 32.0124Dir 


510.407V - 4331.01(T-T, ) 


ent 








TABLE V 
Area of Station Ship E 
35N 48w 


June December 


Number of 
Observations 78 103 
—_" | 


Standard Error 0.4446 G60 1a" 


Mest Sienificant 
eee cies u,V,RH,E Ty oV se 2F 


° 
* 










Regression Equations 
dume: Vis = 865.3962 — 0707962u +001 jeu 

+ 281 {IAN + 00707 5 Cau - 1.68866e, - 0.071146 rem 

- 0.08467E - 0.046262RH- 3.10378WVT + 0.01130Dir 


TOMO Oa tv — 0.13120(T-T,) 


December: Vis = 90.5714 - 0.01468u + 0.04646v 
+ eS ro + 0.23113(T_-T) = 3.8770e - 0.15234 (e,-e.) 
+ 0.20835E - 0.15799RH + 12.1009WVT + 0.00659Dir 


-14194V - 0.64965(T-T,) 


Ce 


TABLE VI 
Area of Station Ship I 
59N 19W 


a 
Number of 81 100 
Observations | 

—_ 

Standard Error 0.9712 
Most Significant V,v,WVT,u 
Variables 













Regression Equations 






Jumes. Vic = 109,024 4 0.07426un- 70) Osome., 







~30439T - 103e9T . - 0.02041e + COO Sale F 


.12608RH + 18.9350WVT - 0.01454Dir - 0.22802V 






CT) 





December: Vis = 102.565 - 0.04026u - 0.37575T 










6293958, + 0.05456E - 0.13815RH + 39.2677WVT 


.O2940Dir - 0.08429V - 0.14386(T-T,) 








TABLE VII 
Area of Station Ship J 


52.5N 20W 
Number of 71 100 
Observations ; , 


0.5864 OR Soule 


aan 
V 


Most Significant 


pve 
Variables 





Regression Equations 


ime: Vis = lol. 007 — 0.242220 — 0. iseee, 
0.797587. ~ SOAS SRM, - 1.17 ae + 3.0717 oer 
eben + 36es4osvevT = 0.15858V - 1 fn Oeeia) 


December: Vis = -— 5840.40 - 576.669u - 86.5089v 


Dat ieoer. + 6165.60e, | + 895.029E + 199165). 9 ae 


18042 .2WVT - 747.738V 
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TABLE VIII 
Area of Station Ship K 
4ON 16W 


—- »..% June December 


Number of 77 Ql 
Observations 


Regression Equations 


dune: Vis = 144747 - 0.15999u + 0.130700 
- 2 O20 Se, - 2.49072e, at 1 eee +; 1.60880(e,-e.) 


“Orde oRnn + O.00949Dar= + 0.25229V = 3.84869 (T-T, ) 


December: Vis = 132.597 - 0.03920u - 0.02028v 


+0.59468(T_-T) + ie Seem = 0.792e7e. - 0. 4233 eran 


- 0.46488RH + 0.01406Dir - 0.01764V - 1.11553(T-T,) 





2) 





TABLE IX 
Area of Station Ship M 
66N O2E 


| ee 
Number of 38 24 
Observations 


“_ 


Mest Stenificant 
Variables 


Regression Equation 


June: Vis = 112.891 + 0.14855u + 0.12889v ~ 


2.18478e, - 2.47101le. - 0.11022(e, -e_) +50. 22238 


OmhoemoRH — 23.50/41 Vev ie - 0.00220Dir - 0.01236V 
0.43712T, = 2.39449 (T-T,) 


December: Vis = 138.993 + 0.25344u + 0.18560v 
8.87585e. = 1.32500(e,-e.) + 0.67601E + 0.17117RH 
Flex0i50'3 ee - 0.00608Dir - 0.24460V + 4.390007, 


4.98978(T-T,). 
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TABLE X 


Norwegian Sea Area 


Number of 26 lh 239 
Observations 


Standard Error 0.9720 


Most Significant 
Variables 













Regression Equations 





June: Vis = 108.569 + 0.01895u - 0.07545v 
- Ono23280T -T) _ 0.42043e, ae 0.27283e. + 0.01940E 
- 0.09583RH - CMO07 Ciba — OSS oey + 0.052531, 


- 0.65299(T-T, ) 


December: Vis = 106.656 + 0.02474u - 0.02059v 
- 0.388577 - 0.04926T_ - 0.3488le, + 1.038108, 
- 0.89120(e,-e,) - O.12991E - O.14243RH - 0.03066Dir 


.03825V - 0.80111(T-T, ) 
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TABLE 
Grand Banks Area 


Number of 501 
Observations ; 


(ees 
Standard Error 10,069 meters 


Most Significant 
Variables 


& 
;4 
WI 
\O 
—~] 
© 
4 
—~] 
OV 
\O 





Reeress LOM aMalysis 


wane ce Vis = 73957.3°4 17 46620 — 203.2230 
654.609T. - 1284 .42e, = (SES oS + 248 .406e, 
681.950E - 620.758RH + 70.4302Dir - 306.076V 
363.880T, = 3138.57 (T-T,,) 


December: Vis = 50109.7 + 157.05lu - 10.1304v 


ANB) L IU 3 Ti 1398.67(T_-T) 51° 19e.729e. - 179.374E 


193.697RH - 62066.2WVT - 509.609V + 1582.03T, 
NOT EN ee, 
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June 
1. Dewpoint spread (T,; - T 
p p ( Fi ry 


2. Relative Humidity (RH) 


3. Windspeed (V), 


December 
1. Windspeed (V) 
2. Relative Humidity (RH) 


Sc Wil 


Forecast visibilities were penerated USINge Che yore 
swmactonm model parameters as input. As expected, results 
eecarmmed using the full set of input parameters proved to 
be much more meaningful than results obtained without tem- 
Peravure, vapor pressure, and related parameters. Scatter— 
diagrams Viememeonstructed wait the forecast visibiliprcce 
enese indieate relatively little skill. The diagrams 
Pevealed the greatest errors to lie in the low-visibility 


range. 


ao 





V. CONCLUSIONS 


The primary parameters determining visibility at sea 
are wind speed and the relative moisture content of the air. 
Mes verifies G.1. Taylor's study of fog in these tau, 
of the Grand Banks in 1917 and the results of Schramm in 
1966 {10 ,1]. 

Based on the differences in regression equations at 
Bierce rent Station Ships at different periods of time ja 
Single regression equation would not be satisfactory over 
the entire North Atlantic. Further study based on Station 
oieeeecaaca only is necessary to determine the time, period: 
paemercal exvent over which a regression equation could ge 
Semsrtdered valid. 

Mie lack of weight assigned to the W-vt. parameter 
Perests that the gradient of the sea surface temperature 
Was poorly defined in some areas. Undoubtedly important 
Small scale features of the sea surface temperature field 
are smoothed out in the 63 x 63 grid size output of the 
Fleet Numerical Weather Central sea surface temperature 
analysis. This is in agreement with Laevastu's 
ewservations [1] . 

Correlation coefficients were Higher in the cases 
where the majority of the observations were those from 
Station Ship data. A comparison of Grand Banks or Norwegian 


sea correlation coefficients with any Station Shipvaved 
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verifies this fact. Further, the large percenvareuomm- 
plained variance (R*) foreseaeen Ship A in June is attri- 
buted to the fact that the majority of the observations are 
from the Station Ship. It can be concluded that estimates 
Oovisibility at sea by relatively untrained observers are 
mou accurate. This suggests that a scalar analysis of ob-— 
served visibilities would be of limited value except those 
made by trained observers. A visibility analysis based on 
the regression equations may prove to be of some use. 

The results indicate that numerical forecasts of marine 
Permeace parameters such as temperature and humidity have not 
mer reached the accuracy necessary for use in this type of 
Meegiction scheme. 

Given accurate observations of visibility at sea, such 
as those availabie from Station Ship data, reasonably reli- 
able regression equations for visibility based on convention- 
al marine surface observations are possible. Until such 
regression equations are derived and the quality of the 
numerical forecasts of the required parameters is improved, 
numerical-statistical prediction of visibility at sea can 


be used only as rough guidance by an experienced forecaster. 
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